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Abstract. The contributions of different radiation mech- 
anisms to the diffuse 7-ray emission of the galactic disk 
are studied in a broad energy region from 10 4 to 10 14 eV. 
Our analysis shows that at energies between 1 and 100 
MeV the radiation is dominated by the bremsstrahlung of 
relatively low energy, typically less than 1 GeV electrons, 
but with a non-negligible contribution from the inverse 
Compton (IC) scattering of higher energy electrons. Also, 
a significant fraction of the radiation observed at energies 
around 1 MeV could be contributed by mildly relativistic 
positrons annihilating "in flight" with the ambient ther- 
mal electrons. At energies from 100 MeV to 100 GeV the 
7-ray flux is dominated by interactions of cosmic ray pro- 
tons and nuclei with the ambient gas through production 
and subsequent decay of secondary 7r°-mesons. The inter- 
pretation of the GeV 7-ray emission of the inner Galaxy 
as a truly diffuse radiation requires a substantially harder 
spectrum of relativistic protons and nuclei in the interstel- 
lar medium compared with the local cosmic ray spectrum 
measured directly in the solar neighborhood. In the very 
high energy domain, E 1 > 100 GeV, the contribution of 
the IC component of radiation may become comparable 
with, or even could exceed the fluxes of 7r°-decay com- 
ponent if the energy spectrum of electrons injected into 
the interstellar medium extends well beyond 1 TeV. An- 
other signature of multi-TeV electrons is the synchrotron 
radiation which could account for a significant fraction of 
the diffuse hard X-ray flux of the galactic ridge. The fu- 
ture detailed studies of spatial and spectral characteristics 
of 7-ray emission of the galactic disk, especially at very 
high energies around 100 GeV by GLAST, and hopefully 
also at TeV energies by planned ground-based instruments 
should provide important insight into the understanding 
of the sites and mechanisms of acceleration of galactic cos- 
mic rays and the character of their propagation in the 
interstellar magnetic fields. 
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1. Introduction 

Diffuse 7-ray emission of the galactic disk carries unique 
information about the fluxes and the spatial distribution 
of galactic cosmic rays (CRs) . Therefore it is believed that 
the solution of the long-standing problem of the origin of 
galactic CRs essentially depends on the success of obser- 
vational gamma-ray astronomy. Indeed, the separation of 
the radiation components associated with the electronic 
and nucleonic components in a broad energy region from 
1 MeV to 100 TeV would allow determination of the fluxes 
and energy spectra of CRs in different parts of the galactic 
disk, and thus would provide an important insight to the 
character of propagation of CRs in the interstellar medium 
(ISM) . A proper understanding of the latter is a necessary 
condition for accurate estimates of the luminosity of the 
Galaxy in CRs. Our present knowledge about the prop- 
agation of CRs is based on conclusions derived from the 
interpretation of the mass composition and the content of 
the secondary anti-particles (positrons, antiprotons) of the 
locally observed CRs. Although rather effective (see e.g. 
Swordy 1993, Strong et al. 1998), this method requires a 
number of model- dependent assumptions. Moreover, it is 
not yet obvious that the locally observed CRs could be 
taken as undisputed representatives of the whole galactic 
population of relativistic particles. For example, Erlykin 
et al. (1998) recently argued that the fluxes of CRs could 
be dominated by a single or few local sources/accelerators. 
This statement is certainly true at least for the observed 
> 1 TeV electrons which suffer severe synchrotron and IC 
energy losses, and thus could reach us, for any reasonable 
diffusion coefficient, only from the sites no farther than a 
few hundred parsecs (Nishimura et al. 1980, Aharonian et 
al. 1995). 

Therefore the diffuse galactic 7-rays seem to be the 
best carriers of information about the production sites and 
propagation of accelerated charged particles in the galac- 
tic disk (see e.g. Ramana Murthy & Wolfendale 1993). It 
should be noted that the diffuse non-thermal synchrotron 
radiation of the ISM at radio and possibly also at X-ray 
wavelengths provide an additional and complementary in- 
formation, but it concerns only the electronic component 
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of CRs in two extreme energy bands below 1 GeV and 
above 100 TeV, respectively. 

The extraction of the truly diffuse 7-ray emission, i.e. 
the radiation produced by CR electrons, protons and nu- 
clei interacting with the ambient interstellar gas and pho- 
ton fields, is not an easy task because of a non-negligible 
contamination due to weak but numerous unresolved dis- 
crete sources. Before the launch of the Compton Gamma 
Ray Observatory (GRO) the information about the diffuse 
galactic 7-ray background was essentially limited to the 
energy region between 100 MeV and few GeV obtained 
by the SAS-2 (Fichtel et al. 1975) and COS B (Mayer- 
Hasselwander et al. 1982) 7-ray missions. The observations 
and theoretical models based on the results of these satel- 
lites were comprehensively reviewed by Bloemen (1989). 
In brief, these data have revealed a good correlation be- 
tween the high energy 7-ray fluxes and the column density 
of the interstellar hydrogen which was a demonstration of 
the existence of a truly diffuse galactic gamma radiation. 

The observations of the diffuse 7-ray background con- 
ducted in 90's by the OSSE, COMPTEL, and EGRET 
detectors aboard Compton GRO resulted in good quality 
data over five decades in energy of 7-rays (see Hunter et 
al. 1997a, and references therein). These results initiated 
extensive theoretical studies of different 7-ray production 
mechanisms in the ISM (e.g. Bertsch et al. 1993, Gillcr 
et al. 1995, Fathoohi et al. 1995, Gralewicz et al. 1997, 
Mori 1997, Porter & Protheroe 1997, Moskalenko & Strong 
1998, Strong et al. 1998, Pohl & Esposito 1998). 

In this paper we report the results of our study of 
different 7-ray production processes in the ISM. We dis- 
cuss the fluxes of diffuse non-thermal galactic radiation 
of both nucleonic and electronic origin in a very broad 
energy region from hard X-rays to ultra-high energy 7- 
rays. Although at the first glance the problem seems to be 
very complicated and confused because of several compet- 
ing production mechanisms, the existing data of diffuse 
galactic 7-radiation do allow rather definite conclusions 
concerning the relative contributions of different produc- 
tion processes in each specific energy band of 7-rays. A 
separate interest represents the diffuse 7-radiation in the 
very high energy (VHE) domain. In this paper we limit 
our study by the diffuse radiation of the inner part of 
Galaxy at 315° < I < 45°. The inner Galaxy is not only 
the experimentally best studied region in diffuse 7-rays, 
but also, it presents a prime interest because of an en- 
hanced, as currently believed, spatial concentration of CR 
sources in the central part of the Galaxy. We emphasize 
the importance of the multiwavelength approach in solu- 
tion of the problem, and predict a range of 7-ray fluxes 
which could be examined by forthcoming satellite-borne 
and ground based 7-ray detectors. 



2. The spectra of cosmic rays in the inner Galaxy 

We consider a model that assumes diffusive propagation 
of relativistic protons/nuclei and electrons in the Galaxy. 
We will compare our model calculations with observations 
of the diffuse 7-radiation of the galactic disc at |6| < 5°. 
Thus, although the halo of galactic CRs may extend up to 
heights of a few kpc (e.g. see Bcrczinsky et al. 1990, Bloe- 
men et al. 1993), we will need to know the mean spectrum 
of CRs only in a region close to the galactic plane. We 
approximate this region as a disk with a half-thickness 
h ~ 1 kpc, and a surface Su 



2Sg, where Sa = nR, 



G' 

and Rq ~ 15 kpc is the mean radius of the Galaxy. 

An important parameter for calculations of the dif- 
fuse 7-radiation is the mean line-of-sight depth l^ in the 
direction of the inner Galaxy, which basically represents 
the central part of the galactic disk, with some radius 
Rd < Rq, where Rq ~ 8.5 kpc is the distance of the Sun 
from the center of the Galaxy. Because the densities of 
the gas and photon fields in the inner Galaxy are gener- 
ally estimated to be significantly higher than at galacto- 
centric distances R > Rq, the inner galactic disk should 
be responsible for most of the diffuse flux detected at low 
galactic latitudes, neglecting even the effect of a possible 
gradient of CR density towards the center of the Galaxy. 
Thus, the mean value of Id can be reasonably estimated 
as ~ 15 kpc. 

The diffusion equation for the energy distribution / = 
f(r,E,t) of relativistic particles can be written in a gen- 
eral form as (Ginzburg & Syrovatskii 1964): 



dt 



= div r (L>grad r /) - div r (uf) 



_d_ 

~d~E 



{Pf) + A[f\, 



(1) 



where u = u(r) is the fluid velocity of the gas contain- 
ing relativistic particles, and P = P(r, E) = —dE/dt de- 
scribes their total energy losses where we include also the 
adiabatic energy loss term P a( jb = divu_E/3 (e.g. Owens 
& Jokipii 1977, Lerche & Schlickeiser 1982). D = D(r,E) 
is the spatial diffusion coefficient, and A[f] is a functional 
standing for various acceleration terms of relativistic par- 
ticles (i.e. the sources of CRs). 

The integration of Eq.(l) over the volume V — 2hSc 
results in a convenient equation for the total energy dis- 
tribution function of particles N(E,t) — Jfd 3 r in the 
Galactic disk at the heights \z\ < h. The term df /dt leads 
to dN/dt. The volume integral of the two first terms in 
the right hand side of Eq.(l) results in: 

/ [div r (£>grad r /) - div r (u/)] d'V = (2) 

Jv 

(p D (e grad r /) ds - (p (eu)/ds . 



Here e is a unit vector perpendicular to the surface el- 
ement ds directed outward from the disk. These terms 
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describe the diffusive and convective escape of particles 
from the disk through its surface Stot ~ 2 Set- 

The first surface integral in Eq.(2) can be simplified 
if we take into account that the CR density at heights 
z < h ~ 1 kpc, where the diffusion dominates (see below), 
may be significantly higher than at z > h, and approxi- 
mate (egrad/) ~ —n(E,t)/Ax, where n(E,t) = f(E,t) 
is the volume averaged energy distribution function of rel- 
ativistic particles, and Ax = Ax(E) is the characteris- 
tic thickness of the transition layer describing the decline 
(i.e. the gradient) of the density of particles at energy E. 
Taking into account that the total energy distribution of 
particles N(E,t) — 2hSQn(E,t), one finds 



* D (e grad/) = - 

J Stot 



nD2S G 
Ax 



N 

Tdif ' 



(3) 



where Tdif has a meaning of a characteristic diffusive es- 
cape time of relativistic electrons from the Galactic disk: 



1 



D(E) 

T dil (E) hAx(E) 



(4) 



Here D(E) corresponds to the mean diffusion coefficient 
Z? r (r, E) on the surface Sq of the disk. 

The second surface integral in Eq.(2) can be reduced 
to the form —N(E,t)/T conv , which describes a convective 
escape of electrons from the disk through its surface due 
to the galactic wind driven by the pressure of CRs and of 
the thermal gas (e.g. Blocmen et al. 1993, Breitschwerdt 
et al. 1993, Zirakashvili et al 1996): 



ah/u 



(5) 



Here u is the wind velocity on the surface of the Galactic 
disk which could reach ~ 50 km/s (Zirakashvili et al. 1996) 
at the height h = 1 kpc. Thus, the mean convective escape 
time of CRs from the Galactic disk can be estimated as 
Tconv -2x 10 7 ayr 

The parameter a in Eq.(5) is the ratio of the mean 
density n(E) of particles in the disk ( i.e. an average over 
—h < z < h ) to their density f(h, E) at the disk surface, 
z = ±h. Therefore one could expect that a > 1 and gener- 
ally it may be energy-dependent as well, a = a(E). Calcu- 
lations of the spatial and energy distribution f(z, E) of the 
galactic CRs in the framework of the diffusion-convection 
model show (Lerche & Schlickeiser 1980; see also Lerche 
& Schlickeiser 1982, Bloemen et al. 1993) that at eleva- 
tions z <C z c (E), where z c (E) is a characteristic height 
of the diffusion dominated region for particles with en- 
ergy E, the spatial density f(z, E) is almost independent 
of z. Thus, at energies E > £?*, where E* is found from 
the equation z c (E) = h, the parameter a ~ 1. In the ap- 
proximation u(z) = vqz for the galactic wind speed and 
D(E) oc E Sl with <5i < 1 for the diffusion coefficient, the 
height z c = ^2D(E)/vq(1 + <5i/6) (see e.g. Bloemen et 
al. 1993). This is basically the height at which the char- 
acteristic time scale of the diffusive propagation ~ z 2 / D 



equals the convection time scale v$ . The energy E* can 
be estimated of order of a few GeV, taking into account 
that at these energies D ~ 10 28 cm 2 /s (e.g. Berezinsky et 
al. 1990) and that vo ~ 50 km/s kpc (e.g. Zirakashvili et 
al. 1996), which result in z c ~ lkpc. 

In the convection dominated region, z 3> z c (E), the 
spatial density of particles starts to decline as f(z,E) oc 
(z/z c )~ K with k ~ 1.3 — 1.4 (see Lerche & Schlick- 
eiser 1992, Bloemen et al. 1993), thus f(h,E)/f(0,E) oc 
E~ kSi ' 2 . Therefore at energies E <C E. t the parameter a 
should gradually increase, and in principle could be ap- 
proximated as a(E) oc (E/E*)~ , with an exponent A 
much smaller than kSi/2 because at these energies the 
mean particle density n(E), as compared with f(0,E), 
should also decline. However, taking into account that A is 
significantly less than 1, and in order to avoid an introduc- 
tion of an additional model parameter, below we approx- 
imate the convective escape time as energy- independent, 
i.e. with a(E) ~ 1 at all energies. For the CR proton com- 
ponent this approximation is well justified, because for 
E* of order of a few GeV the energy region E <C -E* ef- 
fectively corresponds to sub-rclativistic protons which do 
not represent an interest for this study. Such a simplifica- 
tion appears reasonable also for CR electrons, because at 
energies E ^C 1 GeV the spectral modifications of the elec- 
trons are defined mainly by their Coulomb energy losses 
which take place on time scales significantly shorter than 
the escape losses. 

The volume integral of the third term on the right 
side of Eq.(l) describes the volume-averaged energy losses 
of the electrons with the rate P(E). At last, the volume 
integral of the 4th term in Eq.(l) describes the sources 
of accelerated particles Q(E,t) inside the volume V. The 
final equation for the overall distribution of electrons then 
reads: 



ON 
~dt 



d(PN 



, -- + « 

OE T esc 



Here r esc is the "diffusive 
particles: 



m 



i 



i 



7-difOE) 



(6) 
convective" escape time of 



(7) 



The range of actual energy dependence of r osc is limited 
at high energies because the diffusive escape time cannot 
be less than the light travel time h/c. This obvious re- 
quirement formally follows from the condition that the 
characteristic length-scale Ax(E) for spatial gradients in 
the distribution function f(z, E) cannot be less than the 
mean electron scattering path A sc (-E), otherwise the dif- 
fusion approximation implied in Eq.(l) fails. Taking into 
account that for relativistic particles D ~ A sc c/3, from 
Eq.(4) we find that indeed Tdif(-E) > T m i n ~ 3/i/c. The 
diffusive escape time can be then presented in the form 



Tdif (E) - no (E/E 1 



\-S 



(8) 
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where E\q = 10 GeV. Since Tdif increases for decreasing 
E, t csc given by Eq.(7) becomes energy independent be- 
low some E* when Tdif(P*) = r conv . Neglecting at these 
energies r m i n in Eq.(8), from Eq.(7) then follows that in 
the case of a pure power-law approximation for Tdif (E) as 
in Eq.(8) the overall escape time can be presented in the 
form t csc « r conv /[l + (E/E*) 6 ]. 

Note that generally the power-law index 8 in Eq.(8) 
for the diffusive escape time in the Leaky-box type Eq. (6) 
should be smaller than the index 8\ of the diffusion coef- 
ficient D(E) ex E Sl . Since a faster diffusion of more ener- 
getic particles tends to smooth out the gradients of the dis- 
tribution function f(z, E) more effectively, the character- 
istic length-scale Ax(E) would increase with energy. For a 
power-law approximation Ax(E) oc E 2 the index 82 > 0, 
therefore from Eq.(4) follows that 8 = 8\ — 82 < 8\. This 
consideration may help to qualitatively understand a for- 
mal discrepancy between CR spectral modifications due to 
particle propagation effects predicted in the framework of 
simplified Leaky-box models and more accurate diffusion- 
convection propagation models. Calculations for the latter 
models show (e.g. Lerchc & Schlickciser 1980, Bloemen et 
al. 1993) that at elevations z <C z c (E), where the dif- 
fusive propagation dominates over convection, the initial 
power-law spectrum of injected particles, Q(E) oc E~ T ° , 
steepens by a factor oc E~ dl , whereas at z > z c the in- 
crease of the effectiveness of the convective propagation 
results (in the case of negligible energy losses) in a steep- 
ening to only a half of the power-law exponent of the dif- 
fusion coefficient, f(z,E) oc E~ r °~ Sl / 2 . Meanwhile, the 
(diffusive) escape losses in the Leaky-box models result 
in a single power-law spectrum n(E) oc E~ T °~ . Such a 
difference between the predictions of the CR spectra in 
the framework of diffusion-convection models and Leaky- 
box type models could be qualitatively explained, if we 
take into account that n{E) represents the mean particle 
spectrum integrated over z < h, which is therefore con- 
tributed (in fractions varying with E) by both 'diffusion' 
and 'convection' dominated regions. Thus, the power-law 
index <5 for the escape time should be effectively in the 
region 8\/2 < 8 < 8\. 

Assuming a time-dependent injection function Q(E, t), 
Eq. (6) can be used for determination of the overall energy 
distribution of relativistic particles in a general case of a 
non-stationary source. If the energy losses are independent 
of time, the solution to this equation, in terms of spatial 
density functions n — N/V and q = Q/V, reads: 



n(E,t) 



1 



P{E) J 



P(tt)q({t,h) x 



Here the variable £t corresponds to the energy of a particle 
at an instant t\ < t which has the energy E at the instant 
t, and is determined from the equation 



'-'^kmrr (10) 

For a quasi-stationary injection of electrons into the ISM 
on time-scales exceeding the escape time T esc (E) the en- 
ergy distribution of particles becomes time-independent. 
For calculations of the energy distribution of CR pro- 
tons in the Galaxy we take into account the energy losses 
connected with their inelastic interactions with the ISM 
gas, and the adiabatic losses of particles in a gradually 
accelerating wind. In a simple approximation u{z) oc z, 
used e.g. by Lerche & Schlickeiser (1982) and Bloemen et 
al. (1993), the mean adiabatic energy loss term P a db(-E') 
is found after a simple integration 

rh Idu ., _ , „,„^u{h) 



( 1 811 
2S G / -^f(z,E)Az~N(E)E- 
Jo 3 dz 



3/v 



exp 



Ax 



t x T osc(,Sa:J 



d*i 



(9) 



in the volume of the disk with \z\ < h. This expression 
corresponds to P a db = P/fadb with r a db = ih/u{h) ~ 
6 x 10 7 yr for u(h = 1 kpc) ~ 50km/s. 

For the energy losses of particles due to their interac- 
tions with the gas, which contribute to the overall energy 
loss term in Eq.(6), we should use the volume-averaged gas 
density Sh = L nn(z)Az/h 1 where nu = rifji + 2 x riH2 
is the gas density in terms of 'H-atoms'. In order to esti- 
mate tIh, the HI density distribution by Dickey & Lock- 
man (1990) can be taken: a sum of two Gaussians with 
central densities ~ 0.4 and ~ 0.11 cm -3 and FWHMs of 
210 and 530 pc, respectively, and an exponential with the 
central density 0.064 cm -3 and a scale height ~ 400 pc. 
The molecular gas layer can be approximated by a fur- 
ther Gaussian with the mid-plane density 0.3H2/cm 3 and 
dispersion 70 pc (Bloemen 1987). Such a gas density pro- 
file results in the mean hydrogen density nu ~ 0.15 cm -3 
in the region z < 1 kpc. 

An essential process which defines the spectra of both 
CR protons and electrons in the Galaxy is their energy 
dependent escape. Beside this, for calculations of the en- 
ergy distribution of the electron component of CRs we 
take into account the ionization (Coulomb) losses which 
dominate the overall energy losses of CR electrons in 
the ISM at energies below a few 100 MeV, the adia- 
batic losses, the radiative (synchrotron and inverse Comp- 
ton) losses, and the bremsstrahlung losses. Note that the 
brcmsstrahlung loss term has practically the same energy 
dependence as the adiabatic loss term, Pbrem — Erbrcm 
(neglecting a weak logarithmic increase with energy, e.g. 
see Ginzburg 1979), but for Sh — 0.15cm -3 the cooling 
time Tbrom ~3x 10 7 /njj yr exceeds the adiabatic cooling 
time by a factor ~ 3. Therefore for the formation of the 
energy spectra of CR electrons (but not for the radiation 
flux!) the bremsstrahlung losses are by a factor of 3 less ef- 
fective than the adiabatic losses. At energies above several 
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GeV the total energy losses of the electrons are dominated 
by the radiative energy losses due to synchrotron emission 
in the ISM magnetic field of order of several /iG, and due 
to the IC scattering of the electrons on different diffuse 
target photon fields. 

Besides these processes, which should be taken into 
account for calculations of energy distribution of CR elec- 
trons in the disk, in this paper we discuss also a pos- 
sible contribution to the fluxes of diffuse galactic radi- 
ation caused by annihilation of relativistic positrons in 
'flight' with the ambient thermal electrons. In princi- 
ple, it is possible to include also this process into cal- 
culations of n(E), considering the energy distributions 
for e + and e~ separately, and introducing in the equa- 
tion for n + (E) an additional term which describes the 
disappearance of the positrons due to annihilation on a 
timescale T ann (-E). However, this term does not have a sig- 
nificant impact on the formation of the spectrum n+(E), 
because at any energy the annihilation time T ann (E) = 
(fann(-E)wrTH) -1 > (Trr-pcnn) -1 -3x 10 7 yr, for the mean 
density njj — 0.15 cm -3 , is significantly larger than both 
the Coulomb loss time and the escape times involved. 

For both CR protons and electrons injected into in- 
terstellar medium we assume a stationary source function 
(per unit volume) in a 'standard' power-law form with an 
exponential cutoff at some energy Eq: 

q{E) ex £T r ° exp(-£/£ ) • (11) 

The flux of diffuse radiation with energy E y in a given 
direction is defined by the unit volume emissivity g 7 (r, E 7 ) 
integrated along the line of sight: 



J{E 1 



47T 



-dZ 



q 7 (E 7 )l d 

47T 



(12) 



where q {E 1 ) is the mean emissivity, and Id is the charac- 
teristic line-of-sight depth of the emission region. 

It is convenient to describe the flux of 7-rays produced 
at interactions of CRs with the ISM gas by the emissivity 
per 1 H-atom (see e.g. Bloemen 1989). Then the observed 
intensity of 7-rays linearly depends on the column density 
ATjj along the line of sight. It is worth notice that the 
estimate of the mean hydrogen density n' H — Nn/la, which 
defines the mean emissivities along the line of sight in 
direction close to the galactic plane, may be somewhat 
higher (by a factor about 2) than the mean gas density nn 
which should be used for the calculation of the mean (at 
z < 1 kpc) energy distribution of particles in Eq.(9). Such 
an 'enhancement' of n^ should be allowed, and could be 
understood if one takes into account that at low galactic 
latitudes the radiation fluxes due to CR interactions with 
the gas are not equally contributed by the entire z < 1 kpc 
region, but effectively only by a fraction of this region 
close to the galactic plane, with a thickness of a few 100 
pc, where the spatial concentrations of both the relativistic 
particles and (especially) of the ISM gas are higher than 
their respective mean values averaged over z < 1 kpc. For 



calculations of the IC 7-ray fluxes the same characteristic 
energy densities of the diffuse galactic photon fields as for 
the Eq.(9) should be used, because at elevations z ~ 1 kpc 
these densities are still approximately the same as in the 
Galactic plane (e.g. see Chi & Wolfendale 1990). 

3. Diffuse gamma radiation connected with the 
electronic component of CRs 

There are four principal processes of production of non- 
thermal hard X-rays and 7-rays in the ISM by CR elec- 
trons: inverse Compton (IC) scattering, bremsstrahlung, 
annihilation of positrons, as well as synchrotron radiation 
provided that the electrons are accelerated beyond 100 
TeV. In this section we discuss the first three mechanisms; 
the synchrotron radiation of hard X-rays will be discussed 
in Sec. 5 in the context of the IC radiation of highest en- 
ergy electrons. 

3.1. IC gamma rays 

The calculations of the diffuse IC 7-rays require knowl- 
edge of the low-frequency target photons and of the flux 
of electrons in ISM. 

The photon fields which are important for produc- 
tion of IC 7-rays in the ISM are 2.7 K cosmic microwave 
background radiation (MBR), and the diffuse galactic ra- 
diation contributed by the starlight and dust photons 
with peak intensities around 1 /mi and 100 /im, respec- 
tively. While the density of 2.7 K MBR is universal, with 
Wmbr ~ 0.25eV/cm 3 , the densities of diffuse galactic ra- 
diation fields vary from site to site, and actually are model 
dependent. Detailed calculations of Chi & Wolfendale 
(1991) show that the starlight energy density increases 
from the local value wnir — 0.5eV/cm 3 (Mathis et al. 
1983) up to ~ 2.5eV/cm 3 in the central lkpc region of 
the inner Galaxy. For calculations below we use the mean 
value wnir — 1.5eV/cm 3 . The energy density of the FIR 
produced by dust in the galactic plane is more uncertain, 
and is typically estimated from wfir — 0.05 — 0.1 eV/cm 3 
(e.g. Mathis et al. 1983) to w Fm ~ 0.2-0.3 cV/cm 3 (Chi & 
Wolfendale 1991). For calculations of the IC 7-ray fluxes 
below we adopt wfir — 0.2eV/cm 3 . Fortunately, large 
uncertainties in wfir appear not crucial because at all 
7-ray energies the contribution from IC upscattcring of 
2.7 K target photons significantly exceeds the IC fluxes 
produced on FIR (see below). 

Another source of uncertainties in calculations of the 
fluxes of diffuse 7-rays is the lack of independent informa- 
tion about the flux and the spectrum of galactic electrons 
above several GeV. While radio measurements allow def- 
inite conclusions about the average electron flux below 
a few GeV, at higher energies the electron fluxes in the 
Galaxy are in principle model dependent. The standard 
interpretation of the energy spectrum of CRs usually as- 
sumes a uniform and continuous distribution of sources 
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in the Galaxy both in space and time. Whereas for the 
nuclconic component of CRs this approximation can be 
considered as a reasonable working hypothesis, the valid- 
ity of this assumption for the electrons is questionable at 
least for the high energy part of the measured spectrum 
which extends up to 2 TeV (Taira et al. 1993). Because of 
severe radiative losses, the sources of these electrons could 
not be located well beyond a few 100 pc (Nishimura et al. 
1980, Aharonian et al. 1995), and therefore the measured 
electron spectrum might not be applicable for calculations 
of 7-radiation from the distant parts of the galactic disk. 
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Fig. 1. The fluxes of CR electrons near the Sun from 
sources continuously and uniformly distributed in the 
galactic disk calculated in the framework of the diffu- 
sive propagation model for a diffusion coefficient with 
D(lOGeV) = 10 28 cm/s and a power law index Si = 0.6, 
and a power-law injection spectrum of the electrons with 
r o .o = 2.4. The total flux (solid curve, r$ — pc) is decom- 
posed in order to show the contributions from the sources 
located at distances r > tq for different ro indicated near 
the curves. The hatched region corresponds to the esti- 
mate of the mean flux of low energy electrons derived from 
radio data in the direction of galactic poles by Webber et 
al. (1980). The data points shown correspond to the local 
fluxes of CR electrons measured by different groups (for 
details see Atoyan et al. 1995) 



In Fig. 1 the energy spectrum of CR electrons calcu- 
lated assuming a uniform and continuous distribution of 
the sources in the galactic disk (solid curve) is decom- 
posed to show the contributions from sources located at 
distances r > ro for different ro . It is seen that even 
at energies ~ lOGcV the total flux of the observed elec- 
trons is dominated by particles accelerated and injected 
into ISM at distances r < 1 kpc from the Sun. At TeV en- 
ergies the sources beyond 500 pc contribute only ~ 10% 
of the total electron flux. Since for these relatively small 
spatial scales the assumption of continuous distribution 
(both in space and in time) of CR sources may not be 



well justified, the spectrum and the flux of high energy 
electrons at TeV and higher energies may show significant 
variations in different sites of the galactic disk (Atoyan 
et al. 1995). In particular, one could expect a significant 
enhancement of the electron fluxes in the central region of 
the Galaxy due to presumably higher concentration of cos- 
mic ray sources there. Therefore we may allow deviations 
of the predicted electron distribution in the inner Galaxy 
from the observed fluxes, perhaps except for the region 
below few GeV where the radio observations, provide in- 
formation about the average spectrum of galactic electrons 
along the line of sight (see also Porter & Protheroe 1997, 
Pohl & Esposito 1998, Strong et al. 1998). It should be 
noted, however, that because of a significant absorption of 
radio fluxes in the interstellar medium, a direct informa- 
tion about the spectra of CR electrons in the inner Galaxy 
is not actually available. The fluxes of radio electrons in 
the Galaxy are generally deduced using the observations 
from the directions of the galactic poles or the anticenter 
(e.g. see Webber et al. 1980) under an assumption of a 
homogeneous distribution of CRs in the galactic disk. 

In Fig. 2 we show the average spectrum of electrons in 
the inner part of the Galaxy calculated in the framework 
of the model described in section 2, assuming a power- 
law index for the electron injection spectrum r o .o = 2.15, 
and normalizing the energy density w c = J En(E) dE to 
0.05 eV/cm 3 . For the parameters used in calculations, this 
normalization requires an acceleration rate of electrons 



L c ~ 1.6 x 1(T' erg/kpc 3 s . 



(13) 



This implies that for the inner part of the Galactic 
disk, with R < 8.5 kpc and a half thickness 1 kpc, the 
overall acceleration power should be about 7 x 10 39 erg/s. 
In the energy region between 100 MeV and 1 GeV the 
energy losses of electrons are dominated by the adiabatic 
losses and the brcmsstrahlung, with dE/dt oc E, which do 
not change the original (acceleration) spectrum of elec- 
trons. Indeed, it is seen from Fig. 2 that in this energy 
region, which is responsible for synchrotron radio emis- 
sion, the electron spectrum remains rather close to the 
injection spectrum. Therefore the spectral index of the ob- 
served synchrotron radio emission contains almost model- 
independent information about the injection spectrum of 
electrons. At energies below 100 MeV the electron spec- 
trum suffers significant deformation (flattening) because 
of ionization losses, and at energies above 1 GeV the en- 
ergy distribution of the electrons steepens because of a 
combination of escape and radiative (synchrotron and in- 
verse Compton) losses. 

The calculated fluxes of diffuse 7-rays produced by 
electrons are shown in Fig. 3. The spectrum of IC 7- 
rays below the highest energy observed by EGRET, E < 
30 GeV, is not very sensitive to the exact value of the 
cutoff energy Eg in the injection spectrum of electrons, 
provided that Eq exceeds 10 TeV. For the energy density 



F.A. Aharonian and A.M. Atoyan: Diffuse galactic gamma-ray emission 



1000 




1 2 

Log(E/GeV) 

Fig. 2. The mean flux of electrons in the central region 
of the Galaxy calculated assuming a stationary injection 
spectrum of electrons with r o .o = 2.15 and Eq = 100 TeV, 
and the following model parameters for the escape time 
in Eq.(7): r conv = 2 x 10 7 yr, t w = 10 7 yr, 5 = 0.6, 
7"min = 3 x 10 3 yr. For parameters of the ISM we have as- 
sumed: B = 6/iG, nn = 0.15 H — atom/cm 3 (see text), 
and wfir = 0.2eV/cm 3 , wnir = 1.5eV/cm 3 . The in- 
jection rate of electrons is normalized so that the en- 
ergy density in the resulting spectrum of CR electrons 
is w c — 0.05eV/cm 3 . The hatched region is the same as 
in Fig. 1 

of the diffuse interstellar NIR/optical radiation we have 
assumed wnir = 1.5eV/cm 3 . For this value of wnir the 
IC radiation component produced on the galactic starlight 
photons (dot-dashed line) somewhat exceeds in the energy 
region 10 MeV - 30 GeV the IC flux produced on 2.7 K 
MBR (dashed line). The 'FIR' component of IC radiation 
(3-dot-dashed line) calculated for wfir = 0.2eV/cm 3 at 
any 7-ray energy contributes less than 25% of the total IC 
flux. 

In its turn, the overall IC 7-ray flux can account, for 
the chosen infrared photon field densities, only for < 20 % 
of the observed 7-ray fluxes both at MeV ("COMPTEL") 
and GeV ("EGRET") energies (see Fig. 3). For the same 
average electron fluxes shown in Fig. 2, the fluxes of IC 
radiation could be increased assuming formally a larger 
depth Zd of the emission region. However, the value of the 
mean ^d = 15kpc assumed in Fig. 3 is already large, and 
hardly it could be significantly increased further. Another 
way to increase the flux of IC 7-rays is possible if we as- 
sume that the energy density of the electrons in the in- 
ner Galaxy is significantly larger than w c — 0.05eV/cm 3 . 
However, for a fixed gas column density Nh this would 
automatically increase also the flux of the brcmsstrahlung 
7-rays, resulting in an overproduction of diffuse radiation 
in the 10 — 30 MeV region (see below). 

The possibilities to increase the flux of IC 7-rays to 
a level significantly higher than in Fig. 3 are essentially 
limited also by the radio observations. The density of the 
Galactic electrons in the energy range 70 MeV to 1.2 GeV 




Log(E/eV) 

Fig. 3. The flux of diffuse 7-rays produced by the CR 
electrons due to different radiation processes in the in- 
ner Galaxy. For calculations we assume the characteristic 
line-of-sight depth of the emission region Id — 15kpc and 
the gas column density iVn = 2 x 10 22 cm~ 3 . The open 
dots show the bremsstrahlung flux, and the open trian- 
gles show the overall flux of the IC radiation due to dif- 
ferent target photons: 2.7 K MBR (thin dashed line), dif- 
fuse NIR/optical radiation assuming wnir = 1.5eV/cm 3 
(dot-dashed line), diffuse FIR radiation with uipiR = 
0.2eV/cm 3 (3-dot-dashed line). The heavy dotted line 
shows the flux of 7-rays due to annihilation of relativistic 
positrons in flight in the case of a high charge composi- 
tion C+ = e + /(e + + e~) = 0.5 for electrons with energies 
1 — 100 MeV. The sum of the bremsstrahlung and IC fluxes 
is shown by solid line. The heavy dashed line corresponds 
to the overall 7-ray flux including also the annihilation 
radiation. The data points show the mean flux of diffuse 
high energy 7-rays observed by EGRET (Hunter et al. 
1997b), and the hatched region shows the range of aver- 
age diffuse 7-ray fluxes detected by COMPTEL (Strong 
et al. 1997; Hunter et al. 1997a ) from the direction of the 
inner Galaxy at low galactic latitudes. 



shown in Fig. 2 is derived by Webber et al. (1980) from 
radio observations at low frequencies in the galactic pole 
directions assuming the average magnetic field B ~ 6 /zG. 
Thus, below few GeV the spectral index of electrons is well 
fixed, T c = 2 a r — 1 = 2.14 ± 0.06 (for the photon spectral 
index of the observed radio emission a r — 1.57 ± 0.03), 
but their absolute flux depends on the magnetic field. 

The mean energy of the IC 7-rays produced by an elec- 
tron with energy E e on target photons with an energy eo 
is E\c = (4/3) (E c /m c c 2 ) 2 e ■ Therefore IC radiation of 
1 GeV 'radio' electrons on the IR/optical photons with 
eo ~ 1 — 2 eV corresponds to energies E ~ 5 — 10 MeV. 
The expected energy flux, F(E) = E 2 J(E) 1 of these IC 
7-rays can be estimated analytically: 
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for r c ~ 2.15. The comparison of Eq.(14) with the results 
of numerical calculations (the dot-dashed line in Fig. 3) 
shows a reasonable accuracy of this convenient analyt- 
ical expression. Normalizing the flux of GeV electrons 
to the radio flux , F v — vS v ~ 1CP 10 erg/cm 2 sster at 
v = 10 MHz, we find a direct relation between the IC 7-ray 
fluxes produced on NIR/ optical photons by GeV electrons 
and the un-absorbed radio fluxes to be expected from the 
inner Galaxy: 
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For a given radio flux in the galactic plane, a decrease 
of the magnetic field by a factor of two would lead to an 
increase of the electron flux shown in Fig. 2 by a factor 
of 2 Qr ~ 3, and correspondingly to the increase of the IC 
7-ray fluxes by the same factor. 

Unfortunately, at low frequencies the absorption of ra- 
dio fluxes from the inner Galaxy at low galactic latitudes 
is very significant, so the diffuse radio flux Fiomhz is un- 
certain. Eq.(15) predicts that this analytical estimate of 
the IC flux (for E < 10 MeV) would be in agreement with 
the results of numerical calculations shown by the dot- 
dashed line in Fig. 3 for the un-absorbed radio flux from 
the inner Galaxy Fiomhz ~4x 10~ 10 erg/cm 2 sster, i.e. 
by a factor 5 larger than the Galactic radio flux observed 
in the polar directions, « 8 x 10~ n erg/cm 2 sster (Webber 
et al.1980). This implies a size of the halo of the Galac- 
tic CRs extending up to heights <~ 3 kpc (to be compared 
with Id ~ 15 kpc used in Fig. 3), which is in agreement 
with the relevant theoretical predictions in the framework 
of the diffusion models (e.g. Bloemen et al. 1993). 

Fig. 3 and Eq. (14) show that any attempt to explain 
the observed 7-ray fluxes at E ~ 1 — 10 MeV, F ^ s > 
10~ 8 ergem -2 s _1 sr _1 by IC radiation either would re- 
quire an energy density of the NIR/ optical radiation at a 
level of > 5eV/cm 3 , which is much larger than it is gen- 
erally accepted for the density of starlight photons in the 
ISM, or would require a very high flux of the radio elec- 
trons in the inner Galaxy. The latter assumption would 
imply, however, very high radio fluxes, exceeding by more 
than one order of magnitude the flux detected from the 
direction of the Galactic poles, unless we assume a (un- 
rcalistically) low magnetic field, B ~ 1/iG or so, in the 
Galactic plane. Moreover, independently of the strength 
of the interstellar magnetic field, such an assumption of 
high flux of E < 1 GeV electrons leads to a simultaneous 
increase and overproduction of the bremsstrahlung flux 
as well, which would then exceed the 7-ray flux observed 
between 1 and 100 MeV. 



3.2. Electron bremsstrahlung 

Since the bremsstrahlung 7-rays below 1 GeV are pro- 
duced by the same electrons which are responsible also 
for the galactic synchrotron radio emission, the differential 
flux J(E) of this radiation in the region from 30 MeV to 
~ 1 GeV should have a characteristic power-law slope with 
an index coinciding with the spectral index of the radio 
electrons r o ~ 2.1 — 2.2. The results of numerical calcula- 
tions are shown in Fig. 3 by open dots. For a power-law 
spectrum of electrons the energy flux of bremsstrahlung 
7-rays can be calculated analytically: 



F hrem (E) ~ 1.3 x 10" 



O.leV/cm 3 10 22 < 
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erg 
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(16) 



In calculations we have assumed a standard composition 
of the interstellar gas (~ 90 % of the molecular and atomic 
hydrogen, and ~ 10% of helium). 

It is worthwhile to compare the bremsstrahlung flux in 
the energy region E ~ 30 MeV with the 'NIR' component 
of IC flux in the region E ~ 10 MeV since both com- 
ponents arc due to the radiation of the same radio elec- 
trons, (although contributed by two different, low-energy 
and high-energy, parts of the power-law distribution of ra- 
dio electrons, respectively). Assuming for the mean pho- 
ton energy of NIR eo — 1 eV, and r o = 2.15 for the radio 
electrons, from Eqs. (14) and (16) we find: 
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(17) 



Comparison of Eq.(15) with the results of numerical 
calculations in Fig. 3 shows a good accuracy of this ana- 
lytical estimate. 

In the region E < 30 MeV the bremsstrahlung flux is 
due to electrons with E c < 70 MeV (i.e. outside the do- 
main of radio emitting electrons) where ionization losses 
result in a significant flattening of the electron spectra 
n(E c ). This results in a drop of i<brem at 10 MeV by a 
factor of 1.5 compared with Fb rom (30MeV). For wnir = 
1.5eV/cm 3 assumed in Fig. 3 the overall IC flux at 10 
MeV is comparably contributed by both NIR and MBR 
target photons. In the case of a higher density of the dif- 
fuse NIR field in the inner Galaxy the overall IC radiation 
at those 7-ray energies will be dominated by the IC up- 
scattering of the starlight photons. Taking all these effects 
into account, one can conclude from Eq.(17) that at en- 
ergies E ~ 10 MeV the bremsstrahlung should dominate 
the overall diffuse emission observed in the direction of 
the Galactic plane, unless one assumes a very high energy 
density of NIR, wnir > 5 eV/cm 3 as adopted by Strong 
et al. (1998). 
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Independently of the density of the NIR in the inner 
Galaxy a conclusion that the contribution of the brems- 
strahlung to the overall flux of the galactic diffuse 7-ray 
background is large, can be derived from the comparison 
of the fluxes produced by the same 70 MeV < E c < 1 GeV 
electrons in the radio and 7-ray regions. At the photon 
energy 30 MeV this results in 



F brcm (30MeV) ~ 28- 
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(18) 



where Fkjmhz is an un-absorbed flux produced in the 
galactic plane. Since the latter cannot be less (and, pre- 
sumably, is even several times higher) than the flux 
10 -10 erg cm -2 s _1 sr _1 detected from the direction of 
the Galactic poles (Webber et al. 1980), the contribu- 
tion of the bremsstrahlung to the overall diffuse flux of 
E ~ 10 MeV 7-rays should be significant, unless one would 
assume either a very high magnetic field (B ^> 10 /iG) or 
a very low gas column density (Ah <C 10 22 cm~ 2 ) in the 
direction of the inner Galaxy. 

3.3. Annihilation of CR positrons in flight 

Since a significant fraction of the CR electron component 
could be in the form of positrons, a non-negligible contri- 
bution to the diffuse 7-radiation below 10 MeV could be 
due to annihilation of mildly relativistic positrons (Aharo- 
nian & Atoyan 1981a, Aharonian et al. 1983). The differ- 
ential spectrum of the 7-rays produced at the annihilation 
of a fast positron with a Lorentz-factor 7+ = E+/m c c 2 on 
the ambient electrons with density n c is described by a 
simple analytical expression (Aharonian & Atoyan 1981b) 
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where p + — \ n\ — 1 is the dimensionless momentum of 

the annihilating positronFL and the photon energy e = 
E/m c c 2 varies in the limits 



7+ + 1 - P+ < 2e < 7+ + 1 + p+ 



(20) 



For the power-law spectrum of the positrons N + ex 
7^ °, the spectrum of annihilation radiation at high en- 
ergies E 3> m c c 2 has a power-law form 



J ann (£) tx E-^ +1 ^ [ln(2S/m c c 2 ) - 1] 



(21) 



1 Note that in Aharonian & Atoyan (1981b) there is a mis- 
print in Eq.(5), namely the momentum p+ in the denominator 
is missing. 



Thus, the spectrum of annihilation radiation is steeper 
than the spectrum of electrons, in contrast to the spec- 
trum of bremsstrahlung radiation which in the high en- 
ergy limit repeats the spectrum of parent electrons. At 
lower energies the spectrum has a more complicated form 
with a maximum around 1 MeV. The ratio of the fluxes 
J&nn/ Jbrem does not depend on the ambient gas density, 
and for a given ratio C+ = e + /(e + — e~) depends only on 
the spectrum of electrons, being higher for steep electron 
spectra. 

The flux of annihilation radiation calculated for the 
spectrum of electrons shown in Fig. 2, and assuming 50% 
content of positrons, is presented in Fig. 3. It is seen that 
under such an assumption the contribution of the annihila- 
tion radiation at 1 MeV exceeds the fluxes produced by all 
other radiation processes, including the bremsstrahlung. 
Therefore we conclude that depending on the (unknown) 
content of low-energy (< 100 MeV) positrons in CRs, this 
process may result in a significant enhancement of the 
diffuse radiation at MeV energies. Note that at energies 
E c < 1 GeV the fraction of positrons in the local (directly 
measured) component of CR electrons gradually increases, 
reaching C+ > 0.3 (although with large uncertainties) at 
E c ~ 100 MeV (Fanslow et al. 1969). A detailed discus- 
sion of different possibilities which may provide enhanced 
positron flux at low energies is out of the scope of this 
paper. We note only that pulsars could be potential sup- 
pliers of low-energy (e + , e~) pairs into the electronic com- 
ponent of CRs (see e.g. Harding & Ramaty 1987). Also, we 
may speculate that a moderate acceleration of beta-decay 
positrons produced at early stages in SNRs would result 
in such an enhancement. Obviously these possibilities re- 
quire thorough examination, therefore the adopted in this 
paper large content of positrons at low energies should 
be considered as a working hypothesis which helps to ex- 
plain the MeV excess in the galactic diffuse background 
radiation. 

The assumption that the process of annihilation of 
supra-thermal positrons in flight might significantly con- 
tribute to the diffuse low-energy 7-radiation of the in- 
ner Galaxy in principle would imply also a high flux of 
0.511 MeV annihilation line radiation, and a rather broad 
continuum emission at E < 0.511 MeV due to annihila- 
tion of the thermalized positrons through the positron- 
ium channel. OSSE measurements (Purcell et al. 1993) 
have shown that the flux of the diffuse annihilation radia- 
tion detected from the inner Galaxy is dominated (~ 97 %, 
Kinzcr et al. 1996) by the positronium annihilation, which 
consists of two components - a narrow 0.511 annihila- 
tion line and the 7-ray continuum below 0.5 MeV. At 
energies E ~ (0.2 — 0.5) MeV the diffuse 7-ray emis- 
sion of the inner Galaxy is contributed mainly by the 
fluxes of these two components of the annihilation ra- 
diation (not shown in Fig. 3; but see e.g. Hunter et al. 
1997b). The total flux of these photons is at the level 
~ 10 -3 cm _2 s _:l (see Kinzer et al. 1996), which is equiv- 
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alent to ~ 2.4 x 10 2 cm 2 s x ster l for the field of view 
3.8° x 11.4° of the OSSE instrument. 

Our earlier calculations (Aharonian & Atoyan 1981a) 
show that in the approximation of an infinite interstellar 
medium, about 20 % of relativistic positrons would anni- 
hilate in flight on the thermal electrons of the same gas 
medium where they cool due to (predominantly) Coulomb 
and bremsstrahlung energy losses. Thus, in this case the 
total photon flux from the annihilation of the positrons af- 
ter their thcrmalization in the ambient ISM would be by 
a factor of 4 larger than the total photon flux due the an- 
nihilation of relativistic positrons. The integrated flux of 
the annihilation radiation by relativistic electrons shown 
in Fig. 3 (dotted curve), is ~ 6.7 x 10~ 3 cm~ 2 s _1 ster -1 . 
Therefore one could expect that the photon flux associated 
with the thcrmalizcd component of CR positrons would be 
as high as Jo. 511 — 2.7 x 10 _2 cm _2 s _1 ster -1 , provided 
that the positron content in the total flux of CR electrons 
is as high as C+ = 0.5, as assumed in Fig. 3. This flux 
is quite comparable with the annihilation radiation flux 
observed. 

We should note, however, that the ratio of CR 
positrons that annihilate after their thcrmalization to 
the positrons annihilating while remaining still relativis- 
tic may be in fact significantly (by a factor of few) lower 
because this estimate does not take into account the es- 
cape losses of particles from the thin gaseous disk of the 
Galaxy. Meanwhile, the positrons, both relativistic and 
thermalised, escape (in particular, by convection) from the 
disk on timescales comparable with, and even shorter than 
their cooling and annihilation times, which may therefore 
significantly reduce (for the same high C+) the flux of the 
annihilation radiation associated with the thermalized CR 
positrons. This problem needs, however, a separate study 
which is out of the scope of the present paper. 

A possible source of relativistic positrons in the en- 
ergy region below 100 MeV, which presents a prime in- 
terest from the point of view of production of the con- 
tinuum annihilation radiation (by relativistic positrons) 
and the subsequent production of the 0.511 MeV line 
and three-photon positronium continuum (by thermal- 
ized positrons), are interactions of CR protons and nuclei 
with the ambient gas via production and decay of sec- 
ondary 7r + -mesons. The total 7-ray flux above 100 MeV 
from the inner Galaxy, associated with the 7r°-decay com- 
ponent of diffuse radiation, cannot significantly exceed 
~ 10~ 4 ph/cm 2 s sr (see below), therefore it strongly lim- 
its the contribution of these positrons to the observed flux 
of 0.511 MeV line, and consequently also to the 7-ray con- 
tinuum at E < 10 MeV. Obviously a more copious mecha- 
nism for production of positrons in the Galaxy with energy 
< 100 MeV is needed in order to interpret the "MeV" ex- 
cess of the diffuse 7-radiation of the inner Galaxy. Ejection 
of relativistic electron-positron pairs from the pulsar mag- 
netospheres seems an interesting possibility. Apparently, 



this question also needs a separate detailed study, which 
cannot be done in this paper. 

The overall flux of bremsstrahlung and IC diffuse 7- 
rays is shown in Fig. 3 by solid line. It is seen that 
the 7-radiation of the CR electrons is significantly below 
the measured spectrum in the entire energy range from 
100 MeV to 30GeV. Due to the lack of independent in- 
formation on the spectrum of high energy electrons with 
E e ^> 1 GeV in the inner Galaxy, the predictions of diffuse 
7-radiation above several 100 MeV could contain signifi- 
cant uncertainties. Although both at 100 MeV and 30 GeV 
energies the deviation of the calculated fluxes from the 
measurements (by a factor of « 2) should not be overem- 
phasized, the gap by a factor of 5 to 7 around 1 GeV is not 
easy to explain by a reasonable set of model parameters. 
Moreover, the flat shape of the overall flux produced by 
CR electrons cannot explain the GeV bump without vio- 
lation of the fluxes observed at lower energies. Below we 
study the possibility of explanation of this bump by the 
nucleonic component of diffuse radiation connected with 
interactions of CR protons and nuclei with the interstellar 
gas. 

4. Gamma rays of nucleonic origin 

4-.1. Emissivity of it -decay 7 -rays 

Relativistic protons and nuclei produce 7-rays in the in- 
elastic collisions with ambient nucleons due to production 
and decay of 7r°-mesons, pp — > ir° — > 27. This mecha- 
nism has been extensively studied by many authors (e.g. 
Stecker 1979, Dermer 1986, Berezinsky et al. 1993, Mori 
1997). Here we present a simple formalism which allows 
us to calculate with high accuracy the emissivity of 7-rays 
in the case of any broad energy distribution of CRs. 

The emissivity q 1 (E~ f ) of 7-rays due to decay of ir°- 
mesons is directly defined by their emissivity q^{E v ) as 



/■ DC 
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dE v , 



(22) 



where E min 
rest mass. 



E~ + m 2 c /4E 7 , and m^ is the 7r°-meson 



The emissivity of secondary particles from inelastic 
proton-proton interactions can be calculated with high 
accuracy using accelerator measurements of the inclusive 
cross-sections <r(Ei,E p ) for production of a particle i in 
hadronic interactions (see e.g. Gaisser 1990). The emis- 
sivity of 7r°-mesons calculated in the ^-functional approx- 
imation for the cross-section a(E„, E p ) reads 

qA E ir) = cn H / 6(E n - lf w 12 k i n )<7pp(£p)np(.Ep)d.Ep 
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where cr pp (-E p ) is the total cross section of inelastic pp 
collisions, and K v is the mean fraction of the kinetic en- 
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ergy S kin = E p 



D c 2 of the proton transferred to the 



secondary 7r°-meson per collision; n p (E p ) is the energy 
distribution of the protons. 

In a broad region from GeV to TeV energies K^ w 0.17 
which includes also ~ 6% contribution from 77-mcson pro- 
duction (Gaisser 1990). From the threshold at Skin — 
0.3 GeV, Cpp rises rapidly to about 28 — 30 mb at ener- 
gies about Skin < 2 GeV. After that <j pp increases only 
logarithmically. For calculations we approximate 



a pp (E p ) « 30 [0.95 + 0.06 ln(£ kin /l GeV)] mb 



(24) 



for Skin > 1 GeV, and assume a pp = at lower energies. 
More accurate approximation of the cross-section below 1 
GeV (see e.g. Dermer 1986) does not noticeably change 
the fluxes of 7-rays even at very low energies provided 
that the broad power-law spectrum of protons extends 
beyond 10 GeV, and thus the overall flux is contributed 
by protons with energies above few GeV. 

Good accuracy of this simple approach is demon- 
strated in Fig. 4 where we compare the emissivity of ir°- 
decay 7-rays calculated on the base of Eqs.(22)-(24) with 
the results of the recent Monte-Carlo calculations of Mori 
(1997) based on a detailed treatment of the cross-sections 
of secondary pion production at nucleon-nucleon interac- 
tions. The full dots in Fig. 4 correspond to calculations of 
Mori for his "median" proton flux (Eq. 3 in Mori 1997). 
The solid line in Fig. 4 corresponds to our calculations for 
the same "median" proton flux. 

The dashed curve corresponds to our calculations but 
for the local CR proton flux (see e.g. Simpson 1983) in the 
form 



J Q (E P ) = 2.2(S p /lGcV)- 2 ' 75 cm-W^GeV- 1 



, (25) 



which has been used for the detailed 7-ray emissivity cal- 
culations by Dermer (1986). Here also we have a very good 
agreement; the original spectrum from Dermer (1986) is 
not shown in order not to overload the figure with almost 
coinciding curves. 

Two other curves in Fig. 4 correspond to the emis- 
sivities calculated for the power law proton spectra with 
spectral indices T = 2.5 (dot-dashed curve) and T = 2 (3- 
dot-dashed curve), normalized to the same energy density 
of CR protons w p — J n p (E p )E p dE p « 1.2eV/cm 3 de- 
rived from the CR proton flux given by Eq.(25). For CR 
proton spectra with T > 2.4 the emissivities, in terms of 
E'zq^, reach the maximum at E ~ 1 GeV, and then at 
lower energies the spectra sharply decline. Note that we 
do not find any peculiarity in the declining part of the 
spectrum at energies between 100 MeV and 1 GeV nei- 
ther in our nor in the Dermer's (1986) or Mori's (1997) 
spectra, in contrast to the apparent changes of the sign of 
the second derivative in the emissivity spectra presented 
by Pohl & Esposito (1998) and Strong et al. (1998). 
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Fig. 4. The emissivities, per 1 H-atom, of 7r°-decay 7- 
rays calculated using an approximate method given by 
Eqs. (22)— (24) for the spectra of CR protons correspond- 
ing to the 'median' proton flux (solid line) of Mori (1997), 
and the flux given by Eq. (25) (dashed line), as compared 
with the results of detailed emissivity calculations by Mori 
(1997) shown by full dots. The dot-dashed and 3-dot- 
dashed curves correspond to the emissivities calculated 
for the single power-law spectra of protons with indices 
T n = 2.5 and T n = 2. 



4-2. Fitting the GeV bump 

An important feature of the gamma-ray of nucleonic ori- 
gin which allows to fill up the observed 'bump' of diffuse 
radiation at GeV energies in Fig. 3 without violation of 
the fluxes observed at MeV energies is the profound drop 
of the spectrum of this this component below 100 MeV. 

In Fig. 5a we show the fluxes calculated for a single 
power-law spectra of CR protons J(E P ) oc E p p assuming 
different values for T p . All fluxes shown correspond to the 
product of the energy density of protons to column den- 
sity of gas in the inner Galaxy w p A% = 2.5 x 10 22 eV/cm 5 . 
Besides, in order to take into account a contribution due 
to CR nuclei, hereafter the fluxes produced by CR protons 
are multiplied by a constant nuclear enhancement factor 
r/A = 1.5 which takes into account the contribution from 
the nuclei both in CRs and ISM (Dermer 1986), although 
at energies above 100 GeV this factor may gradually in- 
crease (see Mori 1997). 

The dashed curve in Fig. 5a shows that the flux pre- 
dicted for a spectrum similar to the locally observed CR 
protons with the power-law index T p — 2.75, fails to ex- 
plain the diffuse 7-ray flux observed from the inner Galaxy 
at energies above lGeV by a factor of 1.5-2. This deficit 
cannot be removed by an assumption of a larger w p or JVh> 
because the flux predicted at 500 MeV is already equal to 
the observed flux. An assumption of a very hard power 
law index for the CRs in the Galaxy, T p = 2.3 (dot-dashed 
line) , explains the data at few GeV, but over-predicts the 
flux at higher energies. And finally, a moderately steep 
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Log(E/eV) 




Log(E/eV) 

Fig. 5. The fluxes of 7r°-decay 7-rays calculated for CR 
proton energy distributions given in a single power-law 
form (Fig. 5a - top panel), and in the form of Eq.(26) 
(Fig. 5b - bottom panel). The fluxes in Fig. 5a are calcu- 
lated assuming for the product w p Nh = 2.5 x 10 22 eV/cm 5 
and 3 different power-law indices T p : 2.75 (curve 1 - 
dashed), 2.5 (2 - solid), and 2.3 (3 - dot-dashed). The 
fluxes in Fig. 5b are calculated assuming w p iV H = 2.1 x 
10 22 eV/cm 5 , and the same indices T = 2.1 and 6 = 0.65 
in Eq.(26), but three different values for the characteristic 
energy E.y. 3GeV (1 - dashed), 20GeV (2 - solid), and 
100 GeV (3 -dot-dashed). 



spectrum of protons with a power-law index T p ~ 2.5 
could explain the spectral shape of the observed 'GeV 
bump' (solid line). 

There is another possibility to fit the GeV spectrum 
of the observed diffuse radiation by 7r°-decay 7-rays. In 
Fig. 5b we show the fluxes calculated for the spectra of 
CR protons in the form 



n p (E p ) ex E p 



-r 



1 + ' IT 



(26) 



This spectrum corresponds to the power-law index 



of protons IV 



r at energies below some E*, but 



Fp ~ r + 5 at high energies E p 3> E*. An energy distribu- 
tion of protons of this kind is to be expected if the source 
function ( acceleration spectrum) of accelerated protons 
has a single power-law form with an index IV Then the 
energy distribution of the protons, which practically do 
not suffer energy losses (but which nevertheless are for- 



mally taken into account in the numerical calculations), 
can be approximated as n p (E p ) ~ q p (E p ) ■ t csc (E p ). This 
is easily reduced to the form of Eq.(26), if the diffusive 
escape time of CRs from the inner Galactic disk (Eq. 8) 
equals the convective escape time at E p = E*. The re- 
sults of calculations in Fig. 5b show that the GeV bump 
in the spectrum of diffuse radiation can be well explained 
(solid line) by a hard spectrum of accelerated protons with 
To ~ 2.1 if one takes also into account the diffusive escape 
of particles from the galactic disk with a power law in- 
dex S ~ 0.6 resulting in a steepening of the CR spectrum 
above the energy E* ~ 10 - 20 GeV. 

The results shown in Fig. 5 demonstrate that it is pos- 
sible to explain in a natural way the 'GeV bump observed 
in the spectrum of diffuse galactic radiation by 7r°-decay 
7-rays assuming relatively hard spectrum of protons at en- 
ergies below 100 GeV. However, the best fits in Fig. 5 do 
not take into account the fluxes contributed by other 7-ray 
production mechanisms, in particular by the IC radiation 
of the electrons. In fact, the IC fluxes are not negligible 
and should be taken into account for any realistic com- 
bination of parameters characterizing the ISM and CRs. 
It its turn, the expected IC contribution at GeV ener- 
gies is tightly connected with the fluxes at lower energies, 
which are contributed not only by IC 7-rays but also by 
the bremsstrahlung, and possibly also by the annihilation 
components of radiation. Therefore, any attempt to fit the 
observed GeV spectrum of diffuse radiation by 7r°-decay 
7-rays cannot be treated separately, but rather should 
be conducted within the multiwavclcngth approach to the 
problem. 



5. Overall fluxes 

All radiation mechanisms considered in previous sections 
significantly contribute to the overall flux of the broad- 
band diffuse 7-radiation of the Galactic disk. In Fig. 6 
we show that the 7-ray data from ~ 1 MeV to 30 GeV 
can be well explained with a set of quite reasonable pa- 
rameters both for the ISM and for the acceleration and 
propagation of the nuclconic and electronic components 
of the galactic CRs. For calculations in Fig. 6 we have 
assumed a mean line-of-sight depth for the inner Galac- 
tic disk ^d = 15kpc as in Fig. 2, and a column density 
Vh = 1.5 x 10 22 cm -3 . This corresponds to a reasonable 
mean gas density along the line of sight at low galactic 
latitudes about jih = N^/ld = 0.33 cm~ 3 . Calculations 
are done so that the resulting energy distributions of CR 
protons and electrons are normalized to the energy den- 
sities w p — 1 eV/cm 3 and w — 0.075 eV/cm 3 . The latter 
value is larger than the energy density of the local CR 
electrons by a factor about 1.5 or so, depending on poorly 
known flux of the local CR electrons below 1 GeV. Such 
an enhanced energy density of CR electrons is quite pos- 
sible if we take into account that the concentration of CR 
sources presumably increases towards the galactic cetre, 
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and that the electrons suffer significant energy losses. It is 
worth notice that the values of the parameters Nu, l^, w p 
and w e may somewhat vary, but the spectral fits would 
be essentially similar to the one in Fig. 6 if the products 
N-r x w p and «h X w c are kept at the same level as in 
Fig. 6. 




Fig. 6. The fluxes of diffuse radiation produced by both 
electronic and nucleonic components of cosmic rays in the 
inner Galaxy, calculated for a hard power-law source func- 
tions of the electrons with r e o = 2.15, and of the pro- 
tons with r Pi o = 2.1, and the escape time parameters 
S = 0.65, no = 1.4xl0 7 yr, r con v = 2xl0 7 yr. Other model 



parameters are: w c = 0.075 eV/cm 3 



«;,, 



1 eV/cm 3 



N H = 1.5 x 10 22 cm- 2 , ; d = 15kpc, B = 6/zG. Contri- 
butions from 7r°-decay (thin solid line), brcmsstrahlung 
(dashed), inverse Compton (dot-dashed), and positron an- 
nihilation in flight (dotted line, for C+ = 0.5) 7-radiation 
mechanisms are shown. The heavy solid line shows the 
total flux without contribution from the positron anni- 
hilation, and the heavy dashed line takes this flux into 
account. 



In Fig. 6 hard power law source functions in the form of 
Eq.(ll), with r c0 = 2.15 and r p0 = 2.1 for the electrons 
and protons respectively, are assumed. For the CR escape 
times in Eqs. (7) and (8) we have chosen 8 — 0.65, and 
T10 = 1.4 x 10 7 yr and T conv = 2x 10 7 yr. For these param- 
eters the diffusive escape time at the energy E* = 5.8 GeV 
becomes equal to the time of convective escape, so in 
the energy range from a few GeV to ~ 10 GeV the spec- 
trum of CR protons (which do not practically suffer en- 
ergy losses) is described by Eq. (26); it gradually steepens 
from the initial power-law with r p ,o = 2.1 to T p = 2.75. 
A decrease of E*, as compared with the 'best fit' value 
E* = 20 GeV in Fig. 5b, is necessitated by the increasing 
contribution of IC radiation to the observed diffuse flux 
at energies E > 1 GeV. Note that for CR proton spectra 
with T p = r p , + <5 < 2.6 this effect of gradually increasing 
contribution of the IC fluxes at GeV energies makes an in- 
terpretation of the observational data rather problematic, 
whereas without the contribution of the IC component a 



single power-law spectrum of CR protons with T p ~ 2.5 
could explain the 'GeV bump' (see Fig. 5a). 

The heavy solid curve in Fig. 6 shows the total 7-ray 
flux without the contribution from relativistic positron an- 
nihilation. In the energy region below 3 MeV the slope 
of this curve is noticeably flatter than the characteristic 
slope of the fluxes detected by COMPTEL (hatched zone; 
Strong et al. 1997). Meanwhile at MeV energies the flux 
of this annihilation radiation may significantly contribute 
to the overall flux, and it may even exceed the individ- 
ual fluxes of both brcmsstrahlung and IC 7-rays, provided 
that at energies below 100 MeV the positron content in the 
electronic component of the Galactic cosmic rays is signifi- 
cant. The dotted curve in Fig. 6 corresponds to an assump- 
tion of a high value of C + = 0.5, and the heavy dashed line 
shows the overall ("IC+bremsstrahlung+annihilation") 
flux of 7-rays radiated by the electron component of CRs. 
For comparison, we show by stars also the COMPTEL 
data points (Hunter et al. 1997) corrected for the contam- 
ination caused by the positronium annihilation radiation 
observed (Purcell et al. 1993) in the same direction. We 
see that the annihilation of positrons "in flight" , on top 
of the IC and bremsstrahlung fluxes, fits rather well the 
COMPTEL measurements. 

Finally, in the context of principal radiation mecha- 
nisms of the diffuse galactic gamma radiation at MeV 
energies, we should mention also a possible contribu- 
tion of the prompt 7-ray line emission produced by sub- 
relativistic cosmic rays via nuclear de-excitation. The 
cmissivity of the total (unresolved) 7-ray line emission 
in the energy range between several hundred keV and 
several MeV, normalized to the energy density of sub- 
rclativistic CRs w scr = lcV/cm 3 , and calculated for the 
standard cosmic composition (of the ambient matter and 
galactic CRs), is about 2 x 10~ 25 ph/s H-atom (Ramaty 
et al. 1979). This implies that for V H = 1.5 • 10 22 cm~ 2 
and w scr < leV/cm 3 , the energy flux of this component 
of gamma radiation cannot exceed 10 -9 erg/cm 2 sster. 
Therefore, this radiation mechanism could hardly be re- 
sponsible for more than several per cent of the observed 
7-ray flux at MeV energies, unless the energy density of 
sub-relativistic particles in the ISM significantly exceeds 
the "nominal" energy density of CRs in the relativistic 
regime, w scr ^> wq ~ 1 eV/cm 3 . 

Hard X-ray emission from the inner parts of the Galaxy 
has been recently reported by the Ginga and Welcome- 
1 (Yamasaki et al. 1997) , RXTE (Valinia & Marshall 
1998), and OSEE (Kinzcr et al. 1997) teams. In Fig. 6 we 
show the fluxes of diffuse hard X-rays reported by OSSE. 
It is seen that the Bremsstrahlung + IC emission of the 
'radio' electrons may contribute only (10 — 20) % to the 
galactic hard X-ray background. In order to explain the 
bulk of this radiation by the non-thermal bremsstrahlung 
one has to postulate an existence of the interstellar pop- 
ulations of sub-relativistic electrons (e.g. Yamasaki et al. 
1997) and protons (e.g. Boldt 1999). It should be noted 
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that the brcmsstrahlung of sub-relativistic particles, both 
of electrons and protons, is a rather inefficient mecha- 
nism of radiation since due to severe ionization losses 
only ~ 10~ 5 part of the kinetic energy of particles is re- 
leased in the form of non-thermal hard X-rays (Aharo- 
nian et al. 1979, Skibo ct al. 1996). Therefore the "sub- 
relativistic brcmsstrahlung" models require continuous in- 
jection of low-energy electrons and/or protons, e.g. by 
SNRs, into the ISM with uncomfortably large rates, es- 
pecially if one tries to explain the fluxes of the low-energy 
X-rays, W scr ~ 10 43 erg/s (Skibo et al. 1996). Moreover, 
in the case of the proton bremsstrahlung the production 
of X-rays is tightly connected with the prompt 7-ray line 
emission due to the excitation of the nuclei (first of all, Fe, 
C, O, etc.) of the ambient gas by the same protons (Aharo- 
nian et al. 1979). Therefore this allows robust upper limits 
on the flux of the proton brcmsstrahlung X-rays based on 
the observed fluxes of diffuse 7-rays at MeV energies. The 
recent analysis by Pohl (1998), based on the comparison 
of the observed keV and MeV fluxes leads to a conclusion 
that indeed the proton bremsstrahlung alone could not 
be responsible for the bulk of the diffuse Galactic X-ray 
emission. 

An alternative mechanism for an explanation of 
the diffuse X-radiation of the galactic disc - the syn- 
chrotron emission of ultra-relativistic electrons in the 
interstellar magnetic fields - has been suggested by 
Porter & Protheroe (1997). An obvious advantage of this 
mechanism, compared with the bremsstrahlung of sub- 
relativistic particles, is its almost 100% efficiency of trans- 
formation of the kinetic energy of the electrons into the 
X-radiation. On the other hand this mechanism requires, 
for any reasonably ambient magnetic field, an efficient ac- 
celeration of electrons up to energies of 300 — 500 TeV. 

These electrons could hardly be produced by the 
shocks of SNRs because of severe synchrotron losses on 
timescale 



^s 



1.2 • 10 3 (S c /100TeV) _1 (B/10//G)- 



yr. 



(27) 



Comparing this time with the maximum rate of the 
diffusive shock acceleration (Lagage & Cesarsky 1983) in 
the extreme limit of Bohm diffusion, one easily finds an 
estimate of the maximum energy of accelerated electrons 



E 



(max) 



90 



B 



10 nG 



-1/2 



4000 km/s 



TeV. 



(28) 



where v s is the shock speed. Thus, for typical parameters 
of a SNR in a stage close to its Sedov phase, when the 
acceleration of the bulk of relativistic particles takes place 
(e.g. Berezhko & Volk 1997), namely v s < 4000 km/s and 
B > 5 /zG, the characteristic maximum energy of acceler- 
ated electrons cannot significantly exceed 100 TeV. 

More probable sites for acceleration of electrons to en- 
ergies ;§> 100 TeV could be the shocks terminating rela- 
tivistic winds driven by pulsars. For example, it is widely 



believed that the electrons in the Crab Nebula are ac- 
celerated at the termination shock of the pulsar wind to 
energies ~ 10 15 eV (see e.g. Arons 1996). An assumption 
that electrons could be accelerated well beyond 100 TeV at 
the wind termination shocks of much older pulsars could 
then explain the diffuse hard X-ray background radiation 
of the Galaxy as a synchrotron emission of these electrons. 
The life time of E c ~ 300 TeV electrons does not exceed 
several hundred years, therefore they cannot propagate a 
distance more (and probably even much less) than a few 
tens of parsec from their acceleration sites. It means that 
the synchrotron origin of the galactic diffuse X-ray back- 
ground would be possible only in a form of superposition 
of the emission from a number of unresolved weak and 
continuous sources along the line of sight. These condi- 
tions could be satisfied e.g. by 10 5 — 10 6 yr old pulsars, 
the overall number of which in the inner Galaxy could be 
estimated up to ~ 3 x 10 4 if neutron stars are produced 
with a rate of about 1 per 30 yrs. 
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Fig. 7. The diffuse background radiation from the Galac- 
tic plane corresponding to the two-component model for 
the relativistic electrons. The heavy solid line corresponds 
to the fluxes produced by the electrons of the first (main) 
population, with the same model parameters as in Fig. 6, 
but for 5 = 0.7 and C+ = 0.3, and the heavy dashed 
line shows the overall fluxes including the contribution 
from the second ('pulsar wind'- see text) population of 
electrons accelerated to energies E = 250 TeV. The local 
mean magnetic field for the second electron population 
is Bi — 25 /iG. Besides the galactic diffuse background 
7-radiation detected by COMPTEL and EGRET, the 
X-ray backgrounds detected by RXTE (Valinia & Mar- 
shall 1998) and OSSE (Kinzcr et al. 1997), as well as the 
upper flux limits at very high energies (VHE) reported 
by Whipple, HEGRA, and CASA-MIA collaborations are 
also shown. 

In Fig. 7 we show the broad band diffuse radiation 
flux calculated in the framework of the model which as- 
sumes that besides the main population of CRs (presum- 
ably accelerated by the SNR shocks), there is also a sec- 
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ond population of electrons accelerated well beyond 100 
TeV at the pulsar wind termination shocks of old neutron 
stars/pulsars in the inner Galaxy. The parameters for the 
first (main) population of accelerated particles, both the 
electrons and protons, are essentially the same as in Fig. 6 
(except for S = 0.7 and C+ = 0.3 used in Fig. 7). For 
the second electron population we assume an acceleration 
spectrum with T e = 2 and an exponential cutoff energy 
Eq = 250 TeV, but also a turnover of the spectrum at ener- 
gies below 1 TeV. Note that such a turnover at low energies 
is a characteristic feature of the electrons accelerated at 
the pulsar wind termination shocks (see e.g. Arons 1996), 
which however does not affect the overall flux of diffuse 
7-rays below TeV energies (dominated by the radiation 
of the main component of CRs). For the mean magnetic 
field of the ISM in Fig. 7 we assume the same value as 
in Fig. 6, B = B\ = 6/iG. But for the second electron 
component we assume significantly larger ambient field, 
i?2 = 25 fiG, in order to fit the hard X-ray data. Differ- 
ent magnetic fields for the first and second components 
of the synchrotron radiation can be understood if we re- 
member that the multi-TeV electrons cool very rapidly 
therefore they cannot propagate very far from their accel- 
erators. More specifically, the hypothesis that the second 
(ultra-high energy) population of electrons is produced at 
the termination shocks of relativistic pulsar winds (which 
inject not only relativistic electrons, but also magnetic 
fields) could explain also why the magnetic field B^ could 
be much higher than the mean field B\ in the Galaxy. 

The assumption of high magnetic field B2 becomes 
very important for a self-consistent interpretation of the 
'diffuse' flux of hard X-rays detected in the galactic plane 
in terms of synchrotron radiation of E c > 300 TeV elec- 
trons, because otherwise the fluxes of IC radiation pro- 
duced by the same electrons would exceed the upper flux 
limits reported by CASA-MIA collaboration at energies 
E 7 > 100 TeV (Barione et al. 1998). As it is seen in 
Fig. 7, the energy flux of OSSE exceeds the flux upper 
limits of CAS A by a factor of > 300. The ratio of energy 
fluxes of synchrotron to IC radiations produced by an elec- 
tron is about the ratio of the magnetic field to the soft 
photon field energy densities, i.e. (B 2 /87r)/0.25eV/cm 3 = 
10(B/10/iG) 2 , if the IC scattering occurs in the Thom- 
son limit, when the parameter b — 4eoE /(m e c 2 ) 2 <C 1 
(eo ~ 6.5 x 10~ 4 cV is the mean energy of 2.7 K MBR 
photons). Actually, the IC scattering of E c ~ 300 TeV 
electrons occurs in a moderate Klcin-Nishina regime with 
6^3. For these values of parameter b the emissivity of IC 
radiation is reduced (compared with the Thomson limit) 
by a factor 5-10, therefore the ratio of the synchrotron to 
IC fluxes > 300 implies a magnetic field B-i > 20 /xG. 

The acceleration power in the second electron compo- 
nent which is needed for interpretation of the hard X-ray 
background in Fig. 7 is about 6 x 10 36 erg/s per kpc 3 , or 
about L c ~ 1.4 x 10 39 erg/s in the entire inner galactic 
disk with a radius about 8 kpc and a thickness h ~ 1 kpc. 



For the overall number of such pulsars of order 3 x 10 4 , 
this implies a rather modest mean acceleration power per 
1 'old pulsar' of about 5 x 10 34 erg/s, i.e. by four orders of 
magnitude less than the power of the relativistic electron- 
positron wind of the Crab pulsar. Note that the kick ve- 
locities of the pulsars can be of order from a few 100 to 
~ 1000 km/s, so the 10 6 yr old pulsars would be able to 
propagate to distances < (0.3—1) kpc, contributing there- 
fore to the emission at the galactic latitudes up to several 
degree. Note also that even for 10 4 pulsars in the 10° x 90° 
field of view of the inner Galaxy considered here, the mean 
density of such pulsars on the sky corresponds to ~ 10 
per square degree. Therefore in principle for the hard X- 
ray/soft 7-ray detectors, which typically do not have an 
excellent angular resolution, a superposition of the emis- 
sion of large number of such relatively weak sources (with 
sizes presumably > 10 pc ) in the field of view of the de- 
tector could imitate a "diffuse" emission. 



6. Discussion 

The diffuse galactic gamma ray emission carries a unique 
information, a proper understanding of which would even- 
tually result in a quantitative theory of the origin of galac- 
tic cosmic rays. Although the problem is essentially com- 
plicated and confused because of several competing mech- 
anisms of production of diffuse 7-rays, the data obtained 
by the COMPTEL and EGRET detectors aboard Comp- 
ton GRO allow rather definite conclusions concerning the 
relative contributions of different 7-ray production mech- 
anisms in the energy region from 1 MeV to 30 GeV. 

At energies below 100 MeV the diffuse 7-radiation has 
an electronic origin. For a set of reasonable parameters for 
both the ISM and CR electrons, the observed 7-ray fluxes 
from 10 MeV to 100 MeV can be well explained by the 
superposition of the bremsstrahlung and the IC compo- 
nents of radiation of relativistic electrons. At lower ener- 
gies, a non-negligible flux can be contributed by the anni- 
hilation of relativistic positrons with the ambient thermal 
electrons. Moreover, if the content of positrons at ener- 
gies E < 10 MeV is close to e + /(e" + e+) ~ 0.5, the 
7-ray ray flux around 1 MeV would be dominated by the 
non-thermal annihilation radiation. Between several MeV 
and 100 MeV the radiation is dominated by the electron 
bremsstrahlung. This result agrees with previous stud- 
ies (Kniffen & Fichtel 1981, Sacher & Schonfelder 1984, 
Gehrels & Tueller 1993). 

Above 100 MeV the IC process dominates in the pro- 
duction of diffuse 7-rays, but at these energies a new mech- 
anism is needed in order to explain the so called "GeV 
bump" detected by EGRET. We believe that this distinct 
emission feature of the inner Galaxy could be naturally 
explained by 7-rays of nucleonic origin produced at the 
interactions of CR protons and nuclei with the ambient 
interstellar gas. Indeed, the overall observed MeV/GeV 
emission of the inner Galaxy is well fitted by the 7r°-decay 
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7-rays (on top of "IC + brcmsstrahlung + annihilation" 
contribution by CR electrons) provided that the spectrum 
of CR protons in the inner galactic disk at low energies 
is substantially flatter than the spectrum of directly ob- 
served local CRs. In particular, a proton spectrum in the 
form of Eq.(26) with E* ~ lOGcV, T = 2.1 and S ~ 0.6, 
which could be formed due to a reasonable combination 
of diffusive and convective escape time-scales of CRs from 
the inner galactic disk, is able to explain very well the 
observed diffuse galactic 7-ray spectrum (see Fig. 6 and 
Fig. 7) . It is interesting to note that in the total "7r°+IC" 
spectrum the contribution of 7r°-decay 7-ray component 
gradually decreases since the spectrum of CR protons at 



high energies with an index T P 



2.7 is steeper 



than the spectrum of the diffuse gamma radiation with 
a power-law index ~ 2.5 observed at energies above sev- 
eral GeV. However, this decline is compensated by the 
hard IC component of radiation which at energies above 
30 GeV becomes the dominant contributor to the overall 
7-ray flux. For convenience of further discussion, we will 
call this possibility for explanation of the GeV bump in 
the observed 7-ray spectrum the scenario 1. 

Another possibility for explanation of the spectrum of 
7-rays up to 30 GeV could be a single power-law spectrum 
of CR protons with T p ~ 2.5, but then we should assume a 
significant reduction of the IC contribution to the overall 
7-ray flux at such high energies. This possibility, which 
we call the scenario 2, can be realized if the acceleration 
spectrum of CR electrons does not extend to TeV ener- 
gies. In this scenario the 7-ray fluxes below several GeV 
are explained, as in the scenario 1, by superposition of the 
electronic and nucleonic components of radiation. How- 
ever at higher energies the 7-ray flux would be strongly 
dominated by 7r°-decay component of radiation. For the 
acceleration spectrum of protons with To ~ 2.1, the re- 
quired single power-law spectrum of protons in the inner 
Galaxy with T p ~ 2.5 could be formed if S ~ 0.4 and 
E* <C 10 GeV (see Eq. 26). This conclusion would imply 
that the mean time of convective escape of CRs from the 
disk is much larger than the diffusive escape of particles 
at 10 GeV. 

These two scenarios predict essentially different origin 
of 7-rays in the VHE domain. While in the scenario 1 the 
7-ray background of the galactic disk at E > 100 GeV is 
contributed mainly by IC scattering of multi TeV electrons 
on 2.7 K MBR, in the scenario 2 the VHE 7-ray flux is 
dominated by 7r°-decay 7-rays. Therefore the future spec- 
troscopic measurements of the diffuse radiation of the in- 
ner Galaxy at low (|6| < 2°) and high (e.g. 2° < \b\ < 10°) 
galactic latitudes by GLAST at E < 100 GeV, and by 
imaging atmospheric Cherenkov telescope arrays H.E.S.S. 
and CANGAROO-3 (both to be located in the Southern 
Hemisphere) above 100 GeV could provide crucial infor- 
mation about the character of propagation of TeV cosmic 
rays in the galactic disk (Weekes et al. 1997). 



Because all target photon fields for production of IC 
7-rays extend practically with the same intensity (see 
e.g. Chi & Wolfendale 1991) well above the characteristic 
height h ~ 200 pc of the gaseous disk, and since the fluxes 
of very high energy electrons above the galactic plane up 
to several hundred pc could be still significant, the real- 
ization of the scenario 1 can be distinguished from the 
scenario 2 by different spectra of high energy radiation 
detected at different galactic latitudes |b|. In particular, 
the scenario 1 predicts an increase of the relative contri- 
bution of the IC component to the 7-ray flux, therefore we 
should expect some flattening (depending on the rate of 
decline of VHE electron fluxes above the galactic plane) of 
the spectrum of 7-rays at |6| larger than few degree. Note, 
however, that the measurements at different \b\ could not 
distinguish between the scenario 2, which interprets the 
observed Galactic 'excess' GeV radiation in terms of truly 
diffuse emission of (mainly) nucleonic origin from the mod- 
els interpreting this radiation as a superposition of con- 
tributions from unresolved SNRs (Volk 1999, Berezhko & 
Volk 2000), or "active" molecular clouds (Aharonian 1991, 
Aharonian & Atoyan 1996). The level of contamination of 
the truly diffuse radiation of the disk by faint, but nu- 
merous 7-ray sources could be properly estimated only 
by future measurements of the angular distribution of the 
galactic background in scales less than 1°. The GLAST 
with its angular resolution as good as 0.1° and detection 
area by a factor of 10 larger than EGRET (see e.g. Bloom 
1996), is nicely suited for this task. 

The fluxes of 7-rays expected above 1 TeV in both 
scenarios 1 and 2 are below the current flux upper lim- 
its set by the Whipple (Reynolds et al. 1993), HEGRA 
(Schmele 1998) and CASA-MIA groups (see Fig. 8). How- 
ever larger fluxes at such high energies cannot be excluded. 
Indeed, besides the hard 7r°-decay 7-ray component due 
to CR sources which could remain unresolved, we may 
expect also higher IC fluxes of 7-rays radiated by a possi- 
ble second component of ultra high energy electrons with 
a spectrum extending well beyond 100 TeV. The exis- 
tence of such a component of electrons is needed within 
the model which interprets the hard X-radiation of the 
inner galactic disk as a result of synchrotron radiation 
of ultra-relativistic electrons. This hypothesis requires a 
very large (up to -Emax ~ 10 15 eV) maximum energy of 
electrons and large magnetic field B^ > 20 /iG in order to 
fit the observed X-ray fluxes, as well as to avoid an over- 
production of IC 7-rays at 100 TeV. Such large values of 
both -E max arid B2 could be probably explained assuming 
that this energetic component of electrons is produced at 
the wind termination shocks of the ensemble of old pul- 
sars/neutron stars in the galactic disk. The hypothesis of 
the synchrotron origin of hard diffuse X-radiation of the 
galactic disk needs further confirmation based on the de- 
tailed study of spectral and spatial distribution of hard 
X-rays by future satellite missions like ASTRO-E and IN- 
TEGRAL, as well as by a detection of the diffuse 7-ray 
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background of the galactic disk at E ~ 100 TeV at the flux 
level rather close to the current CASA-MIA upper limits. 
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